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■ Abstract 

o 

We calculate the B — > P, 1? — > V (P: light pseudoscalar meson, V: light vector 

P-t. 

O 1 meson) form factors in the large-recoil limit in perturbative QCD approach, including 

d : 

both the vector (axial vector) and tensor operators. In general there are two leading 



components <pB an d 4>b for B meson wave functions. We consider both contributions 
of them. Sudakov effects (k± and threshold resummation) are included to regulate the 
soft end-point singularity. By choosing the hard scale as the maximum virtualities of 
internal particles in the hard 6-quark decay amplitudes, Sudakov factors can effectively 
suppress the long-distance soft contribution. Hard contribution can be dominant in these 
approaches. 
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1 Introduction 



The most difficult task in calculating B meson decay amplitude is to treat the hadronic ma- 
trix element (MiM2\Qi\B), which is generally controlled by soft non-perturbative dynamics 
of QCD. Here Qi is one of the effective low energy transition operators of b quark decays 
[prjl , Mi and M2 are the final state mesons produced in B decays. In the earlier years, this 
hadronic matrix elements of B decays were treated by an approximate method, which is 
called factorization approach |J. In factorization approach the hadronic matrix element of 
four-fermion operator is approximated as a product of the matrix elements of two currents, 
(MiM 2 \Qi\B) ~ (M 1 |j lA1 |0)(M 2 |j2 \B), where j Xft and j% are the two relevant currents which 
can be related to Qi through Qi = ji^j^ • The matrix element of ji M sandwiched between the 
vacuum and meson state M\ directly defines the decay constant of M\. For example, if M\ 
is a pseudoscalar and is V — A current, the relation between the matrix element and the 
decay constant will be (Mi|ji M |0) = ifMiPfi, where /mi and p^ are the decay constant and the 
four- momentum of M 1; respectively. The other matrix element (M2U2 \ B) can be generally 
decomposed into transition form factors of B — > M 2 due to its Lorentz property. The explicit 
definition of B meson transition form factors through such matrix elements can be found in 
section 4 and 5. 

In semi-leptonic decays of B meson, the decay amplitude can be directly related to B meson 
transition form factors without the factorization approximation. For example, for B — > irtvi, 
the decay amplitude can be written in the form: 

A(p B , Pn ) = ^y u6 (%(l - l5 )v e )(n( P7V )\u^b\B(p B )), (1) 

where the form factors F\(q 2 ) and Fo(q 2 ) are defined through the B — > 71 transition matrix 
element (n(p n )\u / ~f fJ 'b\B(pB)) in eq.(0) of section 4. In general the form factors are functions 
of momentum transfer squared q 2 = {p b — P-k) 2 ■ In the region of small recoil, where q 2 is large 
and/or the final particle is heavy enough, the form factors are dominated by soft dynamics, 
which is out of control of perturbative QCD. However, in the large recoil region where q 2 — > 0, 
and when the final particle is light (such as the pion), 5GeV (m B = 5GeV) of energy is 
released. About half of the energy is taken by the light final particle, which suggests that large 
momentum is transferred in this process and the interaction is mainly short- distance. Therefore 
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perturbative QCD can be applied to B to light meson transition form factors in large recoil 
region. 

Before applying perturbative method in this calculation, one must separate soft dynamics 
from hard interactions. This is called factorization in QCD. The factorization theorem has been 
worked out in ref.[|3] based on the earlier works on the applications of perturbative QCD in 
hard exclusive processes ||, where the soft contributions are factorized into wave functions or 
distribution amplitudes of mesons, and the hard part is treated by perturbative QCD. Sudakov 
resummation has been introduced to suppress the long-distance contributions. Recently this 
approach has been well developed and extensively used to analyze B decays || |], [7|, ||, ||, [11], 
[L^, [13| . There is also another direction to prove factorization in the soft-collinear effective theory 
nj, which shows correctly the power counting rules in QCD. In this work we shall calculate 



a set of B — > P and B — > V (P: light pseudoscalar meson, V light vector meson) transition 
form factors in perturbative QCD (PQCD) approach. We use the B wave functions derived 
in the heavy quark limit recently [[HJ, and include Sudakov effects from transverse momentum 
k± and threshold resummation || |l(|. In general there are two Lorentz structures for B wave 
functions. If they are appropriately defined, only one combination gives large contribution, the 
other combination contributes 30%. 

A direct calculation of the one-gluon-exchange diagram for the B meson transition form 
factors suffers singularities from the end-point region of the light cone distribution amplitude 
with a momentum fraction x — > in longitudinal direction. In fact, in the end-point region the 
parton transverse momenta k± are not negligible. After including parton transverse momenta, 
large double logarithmic corrections a s ln 2 kx appear in higher order radiative corrections and 
have to be summed to all orders. In addition to the double logarithm like a s ln 2 k±, there are 
also large logarithms a s ln 2 x which should also be summed to all orders. This is called threshold 



resummation |16| . The relevant Sudakov factors from both k± and threshold resummation can 
cure the end-point singularity which makes the calculation of the hard amplitudes infrared safe. 
We check the perturbative behavior in the calculation of the B meson transition form factors 
and find that with the hard scale appropriately chosen, Sudakov effects can effectively suppress 
the soft dynamics, and the main contribution comes from the perturbative region. 

The content of this paper is as follows. Section 2 is the kinematics and the framework of 
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the PQCD approach used in the calculation of B — > P and B — > V transition form factors. 
Section 3 includes wave functions of _B meson and the light pseudoscalar and vector mesons. 
We give the results of B — > P transition form factors in section 4, and the B — > V transition 
form factors in section 5. Section 6 are the numerical results and discussion. Finally section 7 
is a brief summary. 

2 The Framework 

Here we first give our conventions on kinematics. In light-cone coordinate, the momentum is 
taken in the form k = (^, k±) with k^ — k° ± k 3 and = (A; 1 , k 2 ). The scalar product of 
two arbitrary vectors A and B is A ■ B = A^B* 1 = \(A + B~ + A~ B + ) - A± ■ B±. Our study 
is in the rest frame of B meson. The mass difference of b quark and B meson is negligible in 
the heavy quark limit and we take m& ~ tub in our calculation. The masses of light quarks u, 
d, s and light pseudoscalar mesons are neglected, while the masses of light vector mesons p, u, 
K* are kept in the first order. The momentum of light meson is chosen in the " + " direction. 
Under these conventions, the momentum of B meson is Pb = ^( m B, m>B, 0j_), and in the 
large recoil limit q 2 — > 0, the momentum of the light pseudoscalar meson is Pp = (^|, 0, 0j_). 
For the case of light vector meson, its momentum is Py = ^-(l, 7y, 0j_) with ry defined as 
r v = m v /m B . The longitudinal polarization of the vector meson is El = (^7, — ry, 0±j, its 
transverse polarization St = (0, 0, The light spectator momenta k\ in the B meson and k 2 
in the light meson are parameterized as k\ = (0, ^l^f, ki±) and k 2 = (^2^f , 0, k 2 ±), where k 2 
is dropped because of its smallness (In the meson moving along the 'plus' direction with large 
momentum, the minus component of its parton's momentum k 2 should be very small). We 
also dropped kf because it vanishes in the hard amplitudes, which can be simply shown below. 

The lowest-order diagrams for B to light meson transition form factors are displayed in 
Fig.l. The hard amplitude H are proportional to the propagator of the gluon, i.e., H oc 
1/(^2 — ki) 2 ~ l/(2k 2 ■ ki) ~ l/(k 2 k 1 ). It is obvious that only k{ left in the hard amplitude. 

Factorization is one of the most important part of applying perturbative QCD in hard 
exclusive processes, which separates long-distance dynamics from short- distance dynamics. The 
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Figure 1: Diagrams contributing to the B —>■ P,V form factors, where the cross denoting an 
appropriate gamma matrix. 

factorization formula for B — > P,V transition matrix element can be written as 

dz + d 2 z± dy + d 2 y± 



(P,V(P 2 )\bT tl q'\B(p 1 )) = J d Xl dx 2 d 2 k 1± d 2 k 2± - 



(2tt) 3 (2tt) 3 

x e- ife ^(P,V(P 2 )|g(y)^(0)|0)if^e ifcl - z (0|6(0) p g (J (^)| J B(P 1 )), (2) 

where the matrix elements (P, V(P 2 ) |?(l/)a^(0)|0) and (0\b(0) p q CT (z)\B(Pi)) define the wave 
functions of light pseudoscalar (vector) meson and B meson, which absorb all the soft dynamics. 
H^ a,crp denotes the hard amplitude, which can be treated by perturbative QCD. (3, a, a and 
p are Dirac spinor indices. Both the wave functions and the hard amplitude H are scale 
dependent. This scale is usually called factorization scale. Above this scale, the interaction 
is controlled by hard dynamics, while the interaction below this scale is controlled by soft 
dynamics, which is absorbed into wave functions. The factorization scale is usually taken the 
same as renormalization scale. In practice it is convenient to work in transverse separation 
coordinate space (6-space) rather than the transverse momentum space (&j_-space). So we 
shall make a Fourier transformation / d 2 k±e' l]i± ' h to transform the wave functions and hard 
amplitude into 6-space. 1/6 will appear as a typical factorization scale. As the scale fi > 1/b, 
the interactions are controlled by hard dynamics, and as fi < 1/b soft dynamics dominates 
which is absorbed into wave functions. 

Higher order radiative corrections to wave functions and hard amplitudes generate large 
double logarithms through the overlap of collinear and soft divergence. The infrared divergence 
is absorbed into wave functions. The double logarithms a s ln 2 Pb have been summed to all orders 
to give an exponential Sudakov factor e~ s ^ x ' b,p \ here P is the typical momentum transferred in 



the relevant process, x is the longitudinal momentum fraction carried by the relevant parton. 
The resummation procedure has been analyzed and the result has been given in 0, we do not 
repeat it here. 

In addition to double logarithms a s ln 2 Pb in 6-space (or say fcj_-space equivalently), radiative 
corrections to hard amplitudes also produce large logarithms a s ln 2 x. These double logarithms 
should also be summed to all orders. This threshold resummation leads to || 

S t (x) = 2l+2 ^ /2 + c) [x(l-x)Y, (3) 
Vvrl (1 + c) 

where the parameter c = 0.3. This function is normalized to unity. St{x) vanishes very fast at 
the end-point region x — > and x — > 1. Therefore the factors St(x\) and St (£2) suppress the 
end-point region of meson distribution amplitudes. 



3 The Wave Functions 



In the resummation procedures, the B meson is treated as a heavy-light system. In general, 
the B meson light-cone matrix element can be decomposed as flT], [18| 



iki z 



(2tt) 4 
% 



(0\q a (z)b/3(0)\B(p B )) 



(i)B + rriB)l5 



0a 



i B (k x ) + 



(4) 



/3a 



where n = (1,0, Ot), and v = (0,1, Ot) are the unit vectors pointing to the plus and minus 
directions, respectively. From the above equation, one can see that there are two Lorentz 
structures in the B meson wave function. In general, one should consider both these two 
Lorentz structures in calculations of B meson decays. The light cone distribution amplitudes 



i>3 and 4> B are derived by Kawamura et al. in the heavy quark limit |I5 

Jbx 



0b(x,6) 



/b(A - X 
x/6Ag 



6(A -x)J 



m B bJx(A - x) 



■6(A -x)J 



m B bJ x(A — x) 



(5) 



with Aq = 2A/M B , and A is a free parameter which should be at the order of m B — nib- 



The relation between <p B , <p B and 



'Bi 



arc 



<Pb = <P B > <Pb = (p B ~ 9b- 

The normalization conditions for these two distribution amplitudes are: 



d 4 fci0 B (ki) = 0. 



(6) 



(7) 



From eqs.([S[) and (^), we can see that when x — > 0, 4>b ~h 0; while 0^ ~~ > 0. The behavior 
of 0s with the definition (||) is similar to the one defined in previous PQCD calculations 
IH> S S S i> 0> O' 0- Note that our definition of 4>b, 4>b are different from the previous 
one in the literature [§] 

- 1 d 4 z 



o (2tt) / 



(01^(^)6^(0)15^)) 



v 7 ^ 



($> B + m B )7 5 
S (/t>B + m B )^ 



A 



-^s( k i) + 
^(k x ) + ^#,00 



-M /3a 



Pa 



with 
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B 



(8) 



(9) 



This definition is equivalent to eqs.(§, |]) in the total amplitude. Although the final numerical 
results should be the same, the form factor formulas are simpler using our new definition. 
Another outcome is that our new formula shows explicitly the importance of the leading twist 
contribution <pB, which will be shown later. However, if <p' B and (f)' B are defined as in eqs.([|), 
both of their contributions are equivalently important (see the numerical results in Table 2 of 
Ref . [[TDH ) . It is easy to check that both <p' B and (f)' B here have non-zero endpoint at x — > 0. In 
this case, 4>' B does not correspond to the one defined in the previous PQCD calculations ||, ||, 
where 4>b —>■ 0, at endpoint, when x — > or 1. 

The 7r, K mesons are treated as a light-light system. At the B meson rest frame, the K 
meson (or pion) is moving very fast, one of kf or k{ is zero depending on the definition of the 
z axis. We consider a kaon (or tt meson) moving in the plus direction in this paper. The K 
meson distribution amplitude is defined by 



<K-(P)\s a (z) Uf3 (0)\0> 



dxe 



ixP-. 



75 AM^) + m ~f 5 (f)p(x) - moa^jsPuZv 



(10) 



3a 



For the first and second terms in the above equation, we can easily get the projector of the 
distribution amplitude in the momentum space. However, for the third term we should make 
some effort to transfer it into the momentum space. By using integration by parts for the third 
term, after a few steps, eq.flTOD can be finally changed to be 



<*r(P)|s-(*)« /9 (o)|o> 



V2N~ c Jo 



dxe 1 



75 



K 



x) + m O 750p(x) + m [7 5 (A jb 



x 



da 



(11) 



where ^(x) = and vector n is parallel to the K meson momentum p^- And 

t^ok = fTVjd '(m u + m s ) is a scale characterized by the Chiral perturbation theory. For tt meson, 
the corresponding scale is defined as mo n = m%/(rn u + m^). 

For the light vector meson p, u and K*, we need distinguish their longitudinal polarization 
and transverse polarization. If the K* meson (so as to other vector mesons) is longitudinally 
polarized, we can write its wave function in longitudinal polarization B |2]J 



<K*-(P,e L )\d a (z)up(0)\Q> 

-= J o dxe lxP - z {4 [j> K *<f>K.(x) +m K .<f> K *\ 



x 



+ m K 



(12) 



The second term in the above equation is the leading twist wave function (twist-2), while 
the first and third terms are sub-leading twist (twist-3) wave functions. If the K* meson is 
transversely polarized, its wave function is then 



<K*-(P,e T )\d a (z)u(,(p)\0> = -^=£dxe ixP - z {jl[^ K 



y/2N c Jo 
+im K *e^p a "j^£ v n p v u (i) K *(x)} . 



(f) K *{x) +m K «4> v K «(x) 



(13) 



Here the leading twist wave function for the transversely polarized K* meson is the first term 
which is proportional to 0^-*. 

4 B —> P Form Factors 



The B — > P form factors are defined as following: 



{P(p x )\q llx b\B(p B )) 



mo — mp 

\PB+Pl)n 2 



Fi(q 2 ) + 



'-qMq 2 



(14) 
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where q = Pb — Pi- in order to cancel the poles at q 2 = 0, we must impose the condition 

Fi(0) = F o (0). 

That means in the large recoil limit, we need only calculate one independent form factor for the 
vector current. For the tensor operator, there is also only one independent form factor, which 
is important for the semi-leptonic decay B — ► K£ + £~: 

r 1 2F T (g 2 ) 

{P(pi)\qa ta ,b\B(p B ))=ip 1 q v -q IJ p lu ■ , (15) 

L p J m B + mp 

IF (a 2 ) 

(P( P i)\qa, ul5 b\B(p B )) = enftV-x^- (16) 

Trig + Trip 

In the previous section we have discussed the wave functions of the factorization formula 
in eq.(0). In this section, we will calculate the hard part H. This part involves the current 
operators and the necessary hard gluon connecting the current operator and the spectator quark. 
Since the final results are expressed as integrations of the distribution function variables, we 
will show the whole amplitude for each diagram including wave functions and Sudakov factors. 

There are two types of diagrams contributing to the B —>■ K form factors which are shown 
in Fig.l. The sum of their amplitudes is given as 

F 1 (q 2 = 0) = F (q 2 = 0) 

rl /*oo 

= 8nCpm 2 B / dx\dx 2 / bidb\b 2 db 2 {h e (xi, x 2 , &i, b 2 ) (0b(^i, &i) 
Jo Jo 

(1 + x 2 )4>k(x 2 , b 2 ) +r K (l- 2x 2 ) (<j>x(x 2 , b 2 ) + (f) K (x 2 , 6 2 )) - 0b(xi, 6i) 
4>k( x 2, b 2 ) - r K x 2 (4> K (x 2} b 2 ) + (f) K (x 2} 6 2 )) ) OL s (t\) exp[-Sab(tl)] + (17) 
2r K <f) K (x 2 , b 2 )4> B {xi, bi)a s (tl)h e (x 2 , x u b 2 , bi) exp[-S ab (tl)]j , 

where r K = m 0K /m B = m 2 K /[m B {m s + m d )}; Cp = 4/3 is a color factor. The function h e , the 
scales and the Sudakov factors S a b are displayed at the end of this section. 

For B — » 7r form factors, one need only replace the above K meson distribution amplitudes 
4>k by pion distribution amplitudes 0^ and replace the scale parameter r% by r n = m^jm B = 
m 2 /[m B (m u + md)}, respectively. 

For the tensor operator we get the form factor formulas as 

Fr{q 2 = 0) = 8irCpm 2 B / dx\dx 2 \ bidb]b 2 db 2 {hJxi, x 2 , b\, b 2 ) (0s(xi, b±) 

Jo Jo 
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0i(^2, h) ~ x 2 r K (j>^(x 2 , h) + r K (2 + x 2 )<f) t K (x 2 , b 2 )\ - <Mxi, b x 
4>k(x 2 , b 2 ) - r K $ p K {x 2 , b 2 ) + r K (f) t K (x 2 , b 2 ) ) a s (t\) exp[-S ab (tl)} 
+2r K h e (x 2 , xi, b 2 , 6i)a s (^)0 B (xi, h)(f)^(x 2 , b 2 ) exp[-S ab (t 2 e )]j . 



In the above equations, we have used the assumption that x± « x 2 . Since the light quark 
momentum fraction X\ in B meson is peaked at the small region, while quark momentum 
fraction x 2 of K meson is peaked around 0.5, this is not a bad approximation. The numerical 
results also show that this approximation makes very little difference in the final result. After 
using this approximation, all the diagrams are functions of fcf = X\m B j\[2 of B meson only, 
independent of the variable of kf. 

The function h e , coming from the Fourier transform of the hard amplitude H, is 

h e (x 1 ,x 2 , 61, b 2 ) = K ( v /xiX 2 m B 6 1 ) - b 2 )K (vW"b&i) I {^fx~ 2 m B b 2 ) 

+9(b 2 - bi)K Q {^fx^m B b 2 ) I (v^wb&i)] S t (x 2 ), (19) 

where J is the Bessel function and K , I are modified Bessel functions. 
The Sudakov factors used in the text are defined as 

S ah (t) = s (xim B /y/2, b^j + s (x 2 m B /V2, b 2 j + s ((1 - x 2 )m B /V2, b 2 J 
1 



^ \n(t/A) + ^ ln(t/A) 



(20) 



ln(foiA) -ln(6 2 A). 

where the function s(q,b) are defined in the Appendix A of ref.||. The hard scale Vs in the 
above equations are chosen as the largest scale of the virtualities of internal particles in the 
hard 6-quark decay diagrams, 

t\ = max(- v /^2~mB, I/61, l/b 2 ) , 

t 2 e = max( v /xTm B , I/61, l/b 2 ) . (21) 

5 B — ► V Form Factors 

For the B — > K* form factors, we first define the axial vector current part 
(K*(p 1 )\q^^ 5 b\B(p B )) = i (e* - "^T^m) ( m B + m K *)A 1 {q 2 ) 
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■( ( ^_ s {m 2 B -m? K .) \ A 2 {q 2 ) 



777. _B + WfiT* 



+ « 5 %A(g ), 



(22) 



where e* is the polarization vector of K* meson. To cancel the poles at q 2 = 0, we must have 



2m K *A (0) = (m B + mjf.)^i(°) - (m B - m K *)A 2 (0). 



For the vector current, only one form factor V is defined 

(K\p l )\q 1 ,b\B{p B )) = e, vaP e v *p%p{ 



v * aJ 3 2V(q 2 ) 



(m B + m K *)' 



And for the tensor operators, three form factors are defined: 



(K^p^qa^BiPB)) 



d* drri ( 2 



■rp i 2\ (PB " e *) a 

Vb-V\ 



^P B T 2 {q 2 ) 



(23) 



(24) 



(25) 



(K*( Pl )\qa^j 5 b\B(p B )) = (rfp% - p'ffl 



6T 



Pb-Px 

2\ 



-Uq 2 ) + (^p B -p B nT 2 (q 2 



[e* M rf-J#HW). 



(26) 



Another frequently used set of tensor form factors are defined as below [2^] : 



+ 



e*(m|-m|-.)-(g-e*)(pi+pB)/i ^(g 2 ) 

,2 



Q 



rrin - ml-* 



-(pi 



TJ(^). (27) 



They are useful for the discussion of the flavor changing neutral current decay B — > i^*7 and 
Z? — > K*£ + £~ . The relation between the two set of form factors are 



T[{q 2 ) 
T' 2 {q 2 ) 

n( q 2 ) 



T\(g 2 )+T 2 (g 2 



T 1 (q 2 )+T 2 (q 2 ) + 



T x {q 2 )-T 2 {q 2 



m S — m K* 



(T 2 {q 2 )-T x {q< 



-Ts(q 2 



Pb - Pi 

It is easy to see from the above that at large recoil limit g 2 = 0, T{(0) = T^O). 



(28) 
(29) 

(30) 
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As for B — > p, B — > uo form factors, the definition is similar to the above, just replacing K* 
by p and u respectively. 

Calculating the corresponding amplitude for Fig. 1(a) and (b), we get the formulas for the 
form factors at large recoil as 

A (q 2 = 0) = 8irC F m 2 B dx l dx 2 b 1 db 1 b 2 db 2 x {a s (tl) exp[-S ab (t 1 e )][(j) B (x 1 ,b 1 ) 
((1 + x 2 )(f) K *(x 2} b 2 ) + (1 - 2x 2 )r K * (<j) t K ,(x 2 ,b 2 ) + 0^,(x 2 ,& 2 ))) - 
(4>K*{x 2 ,b 2 ) -x 2 r K * (<p t K *{x 2 ,b 2 ) +0^,(^2,62))) h e (x 1 ,x 2 ,b 1 ,b 2 ) + (31) 
2r K *(j) B (x 1 , bi)(j) s K , (x 2 , b 2 )a s (t 2 e )h e (x 2 , x t , b 2 , b x ) exp[-S ab (t 2 e 



where r^* = rriK*/iTiB- The form factor A Q is the one contributing to the non-leptonic B 
decays B — > PV, where the vector meson is longitudinally polarized. This is shown in the 
above equation (|3l|) that the formula depends only on longitudinal wave functions. On the 
other hand, the form factor Ai contributing to the B — > VV decays depends only on transverse 
wave functions, which is shown below 

Ai(q 2 = 0) = 8nC F m B (m B - m K *) [ dx\dx 2 [ bxdb x b 2 db 2 \h e (xi, x 2 , b±, b 2 ) exp[-S ab (tl)} 

Jo Jo L 

a s (tl) <p B {xiM) (<t> T K *{x 2 ,b 2 ) + (2 + x 2 )r K *(j) v Kt (x 2 ,b 2 ) - r K *x 2 (p a K *(x 2 , b 2 



(p B (x 1 ,b 1 ) U T Kt {x 2 ,b 2 ) +r K *(fr v Kt (x 2 ,b 2 ) - r K *(j) a K ,(x 2 ,b 2 ) 



+ r K * 



(32) 



4>b(xi, bi)[(j) v K ,(x 2 , b 2 ) + (j) a K *{x 2 , b 2 )]a s (t 2 e )h e (x 2 , x h b 2 , bi) exp[-S ab (t 2 e )}} 



The form factor A 2 can be calculated from eq. (|23|) , using the above eqs.(|3~T|, |32"D for A and A\. 
It depends on both transverse and longitudinal wave functions. 

The vector form factor V depending only on transverse wave functions, is expressed as 

V(q 2 = 0) = 8^CFTriB{mB + rriK*) \ dx\dx 2 \ b\db\b 2 db 2 \a s {t l e )h e {xi 1 x 2) bi,b 2 ) 

Jo Jo L 

(V>b(xiA) 4> T K *(x 2 ,b 2 ) + (2 + x 2 )r K *(j) a K *(x 2 ,b 2 ) - r K *x 2 (j) v K *(x 2 ,b 2 ) - 
4> B (x 1 ,b 1 ) 4> T K »(x 2 ,b 2 ) +r K *(J) a Kt (x 2 ,b 2 ) - r K *<\> v K ,{x 2 , b 2 ) ) exp[-S ab (tl)] 
+r K ,h e (x 2 ,x 1 ,b 2 ,b 1 ) [<p v K ,(x 2 ,b 2 ) + (j) a K *{x 2 , b 2 )\ 4> B (xi, h) (33) 
a s (t 2 e ) exp[-S ab (t 2 e )]j . 

As for the tensor form factors, T% and T 2 contributing to the B — > K*^ decay, depend only 
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on transverse wave functions 

T\{q 2 = 0) = 167rC F m| J dxidx 2 b x db x b 2 db 2 ^a s (t 1 e )h e (x 1 , x 2 , h, b 2 ) ^1) 
(1 + x 2 )4> T K ,(x 2 , b 2 ) + 2(1 - x 2 )r K *(p a K «(x 2 , b 2 ) - 2x 2 r K *(p v K *(x 2 ,b 2 ) - 
4> B {xi,bi) <j>K*(x2,b 2 ) + (1 - x 2 )r K *(f> a K*(x 2 ,b 2 ) - (1 + x 2 )r K *(j) v K *(x 2 ,b 2 ) ) 
x exp[-S ab (^)] +r x */i e (x 2 ,xi,6 2 ,&i) [<^*(z 2 , 62) + <j>K*(x 2 , h)] (34) 
x0 B (a;i,6i)a s (tg)exp[-S' a 6(^)]} . 

T 2 (q 2 = 0) = 16-7rC F m 2 B dx x dx 2 \ b 1 db 1 b 2 db 2 a s (t 1 e )h e (x 1 ,x 2 ,b 1 ,b 2 )r K * 

Jo Jo 

[<f>K.(x 2 , 62) -<Pk* (x 2 , 62)] Ub(xiM) -Mxi,bij)exp[-Sab(tl)] . (35) 



While form factor T 3 depends on both longitudinal and transverse wave functions 

^M*? 2 = 0) = 16nC F m 2 B [ dxidx 2 [ bidbxb 2 db 2 \a s (t 1 e )h e (xi,x 2 ,b 1 ,b 2 ) (<p B (x 1 ,bi) 

Jo Jo L 

<Pk* (x 2 , b 2 ) + (2 + x 2 )r K *$ K , (x 2 , b 2 ) - x 2 r K «<\) s K , (x 2 , b 2 ) - 
(j) B (x 1 ,b l ) (f) K *( x 2,b 2 ) +r K *<f>K*(x2,b 2 ) - r K *(f) s Kt (x 2 ,b 2 ) ) 
: exp[-S ab (tl)] +2r K *h e (x 2 ,x 1 ,b 2 ,b 1 )(f) s K *(x 2 ,b 2 ) (36) 
6i)a s (tg) expf-S^tg)]} r^* - 2r^»Ti - T 2 . 



x 

X0b( 



6 Numerical Calculations and Discussions 

In the numerical calculations we use 

A^T 4) = 250MeV, f T = 132MeV, f K = lQOMeV, 

f B = WOMeV, m 0w = lAGeV, m 0K = 1.7GeV, 

M B = 5.2792GeV, f K * = 220MeV, /£. = 180MeV, 

M w = 80AlGeV, f p = 217MeV, fj = lQOMeV, 

f u = 195MeV, fZ = IGOMeV. (37) 

The distribution amplitudes <f> l n (x), (f) l K (x), <f> l p (x) (</> l u (x)) and (f) l K *(x) of the light mesons used 
in the numerical calculation are listed in Appendix A. 
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Figure 2: Contributions to B — ► 7r transition form factors (F B7T (0)) from different ranges of a s /7r, 

(a) with the hard scale chosen as virtualities of internal particles including both quarks and gluons; 

(b) the hard scale chosen as virtualities of only internal gluons. 

Fig.|2|(a) displays the contributions to B — > 7r transition form factor at the large recoil limit 
q 2 — > from different ranges of a s /7t, where the hard scale t is chosen as eq. , i.e., the 
maximum virtuality of both internal quarks and gluons in the hard 6-quark decay diagrams. It 
shows that most of the contribution comes from the range a s / it < 0.3, implying that the average 



scale is around J Aqcd^b- For other B to light meson transition form factor calculations, we 
have very similar results. It is observed that with the hard scale chosen in eq . , PQCD is 
applicable to B — > light meson transition form factors. Recent study shows that PQCD is even 
applicable to B — > form factors ||. However, different perturbative percentage distribution 



over a s /ir was observed in [19[ |22 . We check the reason which causes this difference and find 

that the most important reason is the way of choosing the hard scale t. If the hard scale is 

chosen as the maximum virtuality of only the gluon and other transverse momentum scales, 

i.e., t = max(y i /xiX2mB, l/&i, I/&2), the perturbative percentage distribution will be similar to 

ref . |T9|, 22 1, as shown in Fig.^(b) . Therefore the way of choosing the hard scale is one of the 

important ingredients in the PQCD approach, which deserves more concern. [] Provided that 

the virtuality of the internal quark momentum (longitudinal) must appear as a characteristic 

scale in the hard diagram, in general it should be taken into account. Therefore we think that 
: By numerical check we find that these two different choices of the hard scale only slightly affect the mag- 
nitude of the form factors. For example it can only change the B — > n from factor by a few percent. 
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it is reasonable to choose both the virtualities of internal quarks and gluons as the hard scale. 
Certainly the most powerful proof of this point should be performed under the help of numerical 
calculation of higher order loop corrections. However, such a deeper discussion of this problem 
is beyond the scope of this paper, which shall be left to other attempts. 

Table 1: B meson transition form factors at q 2 = with the hard scale chosen in eq. (|2l|) , and 
the numbers in parentheses are results without the contribution of 4>b ■ 



process 
















B -» 7T 


0.292 ± 0.030 
(0.199) 


0.278 ± 0.028 
(0.189) 












B -> K 


0.321 ± 0.036 
(0.231) 


0.311 ±0.033 
(0.223) 












process 


no) 


4>(o) 


Ai(0) 


A 2 (0) 


Ti(0) 


r 2 (o) 


r 3 (o) 


B^p 


0.318 ± 0.032 


0.366 ± 0.036 


0.25 ±0.02 


0.21 ±0.01 


0.56 ± 0.05 


0.013 ±0.001 


0.06 ±0.01 




(0.226) 


(0.256) 


(0.17) 


(0.14) 


(0.41) 


(0.004) 


(0.05) 


B ^uj 


0.305 ± 0.030 


0.347 ± 0.036 


0.24 ±0.02 


0.20 ±0.02 


0.53 ± 0.05 


0.012 ±0.001 


0.06 ±0.01 




(0.212) 


(0.250) 


(0.16) 


(0.13) 


(0.38) 


(0.003) 


(0.05) 


B -> K* 


0.406 ± 0.042 


0.455 ± 0.047 


0.30 ±0.03 


0.24 ±0.02 


0.69 ± 0.08 


0.007 ±0.001 


0.09 ±0.01 




(0.293) 


(0.336) 


(0.21) 


(0.16) 


(0.51) 


(-0.001) 


(0.07) 



The results of B — > P, V light meson transition form factors are given in Table [I] with the 
hard scale chosen in eq.(^TJ). Compared with previous PQCD calculations on some B — ► P, V 
transition form factors || [| ^| , the current work is different from them mainly on two points: 

1. The B meson wave function used here is what derived from the equation of motion in 
Heavy Quark Effective Theory ||15|| . There is only one free parameter in the functions of 
distribution amplitudes, A. We show the results for A = (700 ± 50)MeV in table [l]. All 
the form factors are sensitive to this parameter, i.e. sensitive to the shape of the B meson 
distribution amplitudes. 

2. Two Lorentz structure terms of B meson wave function, both <pB and <f)B defined in 
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eqs.(B, 0), are taken into account in this work. To see how large the 0# term contributes, 
we give the results without the contribution of 4>b in the parentheses in Table 1. They 
show that the contribution of <f>B is about 30%. The dominant contribution comes from 
(f>B term. This result shows that simply dropping the contribution of 4>b is not a good 
approximation. 

Table 2: Form factors at q 2 = for B — > 7r and B — > p transitions calculated in this work and 
UKQCD . 





F (0) = Fi(0) 


7(0) 


A (0) 


Ax(0) 


A 2 (0) 


UKQCD m 


0.27 ±0.11 


u - oo -0.05 


o 3n +0 06 

u - JU -0.04 


n 9 7+0.05 
u - z/ -0.04 


u - ZD -0.03 


this work 


0.292 


0.318 


0.366 


0.250 


0.210 



We compare some of the results calculated in this work with lattice calculation by UKQCD 



collaboration |25fl in Table 2. It shows that our results are consistent with theirs. 

The B — ► if* form factors are useful for the calculation of flavor changing neutral current 
process B — > K*^ and B — ► K*£ + £~, which have been discussed many times PB[ . We show 
some of them in table |^ for comparison. It is easy to see that our results agree with lattice 



calculations |25| and the results calculated using lattice-constrained dispersion quark model 



7 Summary 

We have calculated B — > P and B — > V transition form factors in the PQCD approach. We not 
only calculate the B to light meson transition form factors defined in vector and axial vector 
currents, but also form factors defined in tensor currents qa^b and qcr^j^b, which can be used 
to study semi-leptonic and radiative B decays induced by magnetic penguin operators qa^l + 
75) bF^p. With the hard scale appropriately chosen, Sudakov effects can effectively suppress 
long-distance dynamics, which makes short-distance contribution dominate. The characteristic 
scale in B to light meson transition processes is around J AQCD m B- 



16 



Table 3: Some form factors at q 2 = for B — > X* transitions calculated in this work and some 
other works. 





T[(o)=mo) 


A)(0) 


A 1 (0) 


quark model [27 


] 


0.155 


0.32 


0.26 


QCD sum rule @ 


0.19 ±0.03 


0.3 ± 0.03 


0.37 ±0.03 


light cone sum rule 28] 


0.18 


0.27 


0.36 


Lattice QCD |25 




lfi+ 02 
u - iD -o.oi 


0.33 


0.29 


dispersion quark model 29] 


0.177 


0.44 


0.33 


this work 


0.175 


0.455 


0.297 



We considered both of the two Lorentz structures of B meson wave functions, and found 
that the contribution of 4>b defined in eqs.(f|, ||) is about 30%. 

Finally we compared our results with Lattice calculation, some quark model and QCD sum 
rule calculations, we found that they are consistent with our results. 
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A Wave Functions of Light Mesons Used in the Numer- 
ical Calculation 

For the light meson wave function, we neglect the b dependence part, which is not important 
in numerical analysis. 

The distribution amplitude (p^ for the twist-2 wave function and the distribution amplitudes 
and $\ of twist-3 wave functions are taken from f20|j , 



<f,A(x) = ^x(l-ar) l±0.44C| /z (t)±0.25C , 4 / "(t) , (38) 



3/ w 
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2V6 
u 

2V6 



/:; r l + 0.43C 2 1/2 (t) + 0.09C 4 1/2 (t) 
1 + 0.55(10s 2 -lOx + 1 



(39) 
(40) 



(41) 



where t = 1 — 2s. The Gegenbauer polynomials are defined by 

Cl /2 (t) = ±(3£ 2 - 1), C\ /2 (t) = §(35£ 4 - 30t 2 + 3), 
Cl ,2 (t) = |(5t 2 - 1), Cf(t) = f (21t 4 - 14t 2 + 1), 

whose coefficients correspond to m 07r = 1.4 GeV. 

We choose the different distribution amplitudes of p meson longitudinal wave function as 



<j) p (x) 



3/ P 



l + 0.18C 2 3/2 (t) 



(3t 2 + 0.3t 2 \bt 2 - 3l + 0.21 [3 - 30t 2 + 35t 4 l } , 



3/J 



1 + 0.76(10x 2 - 10x+ 1)" 



For the transverse p meson we use [21]: 

"l + 0.2C 2 3/2 (£) 



(42) 
(43) 
(44) 



lTi \ 3/J 

rX( 1 — X) 



[X = 



V6' 



V u 

= ^= |^(1 + t 2 ) + 0.24(3t 2 - 1) + 0.12(3 - 30t 2 + 35t 4 )| , 



3 f r 

>°(x) = — -j= t [l + 0.93(10a; 2 - 10s + 1) 



(45) 
(46) 
(47) 



For the u meson, we use the same as the above p meson, except changing the decay constant 
f p with f u . 

We use fix of the K meson twist-2 wave function and 0^ and of the twist-3 wave 
functions from p0| , |21 , |3~0"|[ , 

4>i(x) 
4>Ux) 



3/ 



*x(l -x)\l + 0.51* + 0.3{5i 2 - 1} 



\l + 0.12(3t 2 - 1) - 0.12(3 - 30t 2 + 35t 4 )/8 
2y6 L 

he 
2V6 



l + 0.35(10:r 2 - 10a? +1) 



(48) 
(49) 
(50) 



whose coefficients correspond to m 0K — 1.7 GeV. 
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We choose the light cone distribution amplitudes of K* meson longitudinal wave function 
as El, 



J K' 



X 



X 



fT 

JK* 



x(l — x) 



1 + 0.57t + 0.07C 2 3/2 (t) 

,1/2 



j^j= {o.3£(3t 2 + lOt - 1) + 1.68C 4 1/2 (t) + 0.06t 2 (5t 2 - 3) 
+0.36[1 -2t- 2tln(l -a:)]}, 



-- {3t [l + 0.2* + 0.6(10x 2 - Wx + 1)] - 0.12x(l - x) 

+0.36[l-6ar-21n(l-:r)]}. 

The light cone distribution amplitudes of K* transverse wave function are used as 

3/J 



<Pk*(x) 
<f>Ux) 



X 



— j£x{l-x) 1 + 0.6t + 0.04C 2 3/2 (t) 
v6 L 

0= |^(1 + t 2 + 0.44t 3 ) + 0.4C 2 1/2 (t) + 0.88Q 1/2 (t) 
+0.48[2x + ln(l -x)}}, 

^{341 + 0.19^0.81(10^-10^1)1-1.14,(1-., 
+0.48[l-6x-21n(l-x)]}. 



(51) 



(52) 



(53) 



(54) 



(55) 



(56) 
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